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Abstract 

In this article in a generalization of our previous work, [lj] , we investigate the dynamics of the 
non-local F(R) gravity after casting it into local form. The non-singular bouncing behavior and 
quintom model of dark energy are achieved without involving negative kinetic energy fields. Two 
cosmological tests are performed to constrain the model parameters. In case of phantom crossing 
the distance modulus predicted by the model best-fits the observational data. In comparison with 
the CPL parametrization for drift velocity, the model in some redshift intervals is in good agreement 
with the data. 
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1. INTRODUCTION 



The cosmological observations have provided increasing evidence that our universe is 
undergoing a late-time acceleration expansion j^j. To explain this observation, a variety 
of theories are put forward to introduce a new fluid called dark energy, which possesses a 
negative enough pressure According to the observational data from the Type la Super- 
novae 4] and WMAP satellite j^J, our universe is spatially flat and consist approximately 
of %25 dark matter and %75 dark energy. From the simplest models for dark energy as 
cosmological constant to the recent proposed models such as quintessence 6], phantom 
or combination of these two in a unified model named quintom 8J all try to explain the 
accelerating expansion universe. 

An alternative approach to explain dark energy is the modification of general relativity 
(GR) Here, clS cL generalization of the work in 111, we consider a non-local modified 

gravity, where in addition an arbitrary function of Ricci scalar, F(R), is added in the action, 
with the aim to explain the dark energy responsible for current universe acceleration. For 
cosmic acceleration, the equation of state (EOS) parameter is less that —1/3, while for 
quintessence, Chaplygin gas [l^ and holographic dark energy models 15], 00 always stays 



bigger than —1 and for the phantom models is smaller than —1 [16[. For a dynamical dark 
energy, the EoS parameter evolves from w > -1 to w < -1 (or vice versa) and is called the 
phantom divide crossing. 

In particular, the quintom model of were initially proposed to obtain a model of dark 
energy with the u > — 1 in the past and u < — 1 at present. The model can be viewed as 
dynamical model for dark energy with the feature that its EoS parameter can smoothly cross 
over the cosmological constant barrier u — — 1. To construct a quintum model it is necessary 
to add extra degrees of freedom with un-conventional features to the conventional single field 



theory if we expect to realize a viable model in the framework of gravity theory [18J. Our 
model involves two scalars with one being Quintessence-like and another Phantom-like. The 
model provides a bouncing cosmology which allows us to avoid the problem of the initial 
singularity. 

The paper is arranged as follows: in the next section we establish our model and give 
the exact solution of the scalar-tensor field equations. We also analytically investigate the 
conditions for the EoS parameter in the model to across —1 In section three, by using the 
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field equations from the previous section we obtain numerical solutions for the cosmological 
scale factor, the EoS, Hubble and deceleration parameters. In Section four, we perform 
two cosmological tests to validate our model with the observational data. A summary is 
presented in section five. 



2. THE MODEL 



In this section we consider local formulation of non-local F(R) gravity with the action 
given by fl9|. 

"1 



S 



d XyJ—g 



{R(l + d^d^} + F(R) 



(1) 



where / is a function of scalar field <fr, F is a function of scalar curvature R and ip is a second 
scalar field in the model. Variation of Eq. ([1]) with respect to the metric tensor g^ u gives, 

+ l -{d^d u ct> + d^d^) - (g^n - V M V„)(/ - V) 

- 2( 9flu n - V^ U )F' - 2R IMI/ F' + 9fU/ F, (2) 

where F 1 = dF J^ ■ In a FRW cosmological model, for only time dependent and ip, by 
invariance of the action under changing fields and vanishing variations at the boundary, the 
equations of motion for scalar fields and ip become, 

= 4> + 3H<fi + R, (3) 
= i) + 3HtP - RfW, (4) 

where = d , R = 12H 2 + QH and H is Hubble parameter. Also, the 00 and ii 
components of equation ((2]) have the following form 







.3^(1 + / _ ^) + _ 3H(f^<P — ip) — F 



+ 6(H 2 + H)F' - 36(AH 2 H + HH)F", 
= (2H + 3H 2 )(1 + f-^)-2(H + 3H 2 )F' + ^ 

+ 4 +2H j^ f -^ 2F '^ F - 



(5) 



(6) 



Equations ([5} and (jSJ) can be rewritten as 

3H 2 = {-^0 - 3H(f {1) <p - ip) - F + 6(# 2 + H )F' 

- 36(4H 2 H + HH)F"}/(1 + f -i[>), (7) 
-2H - 3H 2 = {(^ + 2Hj t )(f - j, + 2F>) + l -U 

- 2(H + 3H 2 )F' + F}/(1 + f-i>)- (8) 

In comparison with the standard Friedman equations, one can treat the model as standard 
model with the effect of the gravity modification gives rise to a new energy density and 
pressure for the Friedman equations. So the right hand side of the equations ([7]) and (jHJ) are 

Peff = {^<p-3H(f^<p-ij)-F + 6(H 2 + H)F' 

- 36(4H 2 H + HH)F"}/{1 + f - j>), (9) 
Peff = { l -U+{^ + ^Hj t ){f-^ + 2F')+F 

- 2(H + 3H 2 )F'}/(l + f-^)- (10) 

After some algebraic calculation, we can read the effective energy density and pressure from 
the above equations as, 

Peff = A2 ^4p5 i\M>(A - AF>) - 3HB{A + F') 

- F(A + 2F') -3BF'}, (11) 

Peff = A2 _ X 4F/2 + FB + 2(A + \f')HB 

+ {A + F')B}. (12) 

where A = 1 + / — ijj and B = A + 2F' and p e ff = co e ffPeff- From Eqs. (jSJ) and ( fTTl) we 
obtain, 



H = 2( ^/ 4F , 2) f -C ± ^ - | (^ - 4F^) 2 (2? + 3£F') 1 (13) 

where C = - -0 + 2^)(.4+ F'), X> = -|^0(„4 - 4F') + Fi3 and 8 = / (2) 2 + - 
$ + 2^. Using Eqs. (15]), fTTT]), (112]) and (|I3)) ; one can find, 

H = ^{U-HB + B}- (14) 
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Fig. 1: Plots of the evolution of scale factor a(t) and Hubble parameter H (t) 
at the bounce for F(R) = R 2 , f = fo e fe< ^, /o = —0.5 and 6 = 5. Initial values are 
0(0) = -5, 0(0) = 15, V(0) = 15 , ip(0) = 10, H(0) = 8 and H(0) = 1. 



Also from Eq. (|T 



we have 
H 



1 ( f d 3 d 2 ■ d , 
mi 2 \l*- H ln +H Tt F 



2 (23 - 6/W) L " V^ 3 " ^ 2 ' " ^ 
+ (V> - 4F/ (1) + 3/ (2) 0) + ^(0 + 4if) 



+ 



f 



(3)^2 _ y(2) l H( p + Q H+ 12H' 



+ 2ifU-/ (1) (24#+l)) } 



where £> 7^ 6/' is a constraint on iJ. To explore the possibility of the effective EoS parameter, 
co e ff, crosses —1, the condition, H = —u e ffp)\ u ..^.^.i ^ has to be satisfied. 



3. NUMERICAL RESULTS 



A numerical solution for the corresponding equations for EoS parameter, cu e //, the Hub- 
ble parameter, H, and the scale factor a in our model is presented in this section. In here, 
we choose t = to be the bouncing point. Our solution for H(t), Eq. (fi~3l) . and consequently 
a(t) provides a dynamical universe with contraction for t < 0, bouncing at t — and then 
expansion for t > 0. The above analysis clearly can be seen in the numerical calculation 
given in Fig. 1. It also shows that at the bouncing point t — 0, the scale factor does 
not vanish which means the universe avoid singularity faced in the usual Einstein cosmology. 
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Our solutions for p e ff and p e ff, Eqs. ( TTTT) . (fT2|) shows that the effective EoS parameter 
u e ff approaches negative values less than —1 at about t = and crosses —1 around this 
point as shown in Fig. 2. It also shows that in this model the EoS parameter crosses — 1 



line from u < -1 to w > -1 which is supported by observations 



20). 



In some sense this model is similar to quintom dark energy models that consist two 



quintessence and phantom fields [21] . A detailed examination on the necessary conditions 
requires for a successful bounce shows that during the contracting phase, the scale factor 
a(t) is decreasing, i.e., a < 0, and in the expanding phase we have a > 0. At the bouncing 
point, a — 0, and so around this point a > for a period of time. Equivalently in the 
bouncing cosmology the Hubble parameter H runs across zero from H < to H > and 
H — at the bouncing point. A successful bounce requires that the following condition 
should be satisfied around bouncing point, 

ff = -2^< 1+ ^ >0 ' < 15) 

From Fig. 1 and 2, we see that at t — > 0, u e ff < — 1 and H is positive which satisfies the 
above condition. Also from Fig. 2, around t = 0, for deceleration parameter, q(t), we have 
a large acceleration which followed by a period of deceleration and a second acceleration. 
It should be noted that our model violate the strong energy condition and the null energy 
condition (NEC), (p e // +Peff) < for a short period around the bounce point and maybe 
physically unstable. Moreover, after the bounce the equation-of-state (EoS) parameter is 
able to transit from u e ff < — 1 to u e ff > — 1 and then connect with normal expanding 
history. 

From Fig. 3, it is shown that the numerical solutions of the proposed model can re- 



produce the expected red-shift behaviors of the EoS parameter cu e ff(z) — — 1 at z — 0.29 [22]. 



4. COSMOLOGICAL TESTS 



We now examine our model with the observational data using the following cosmological 



tests [23|-|24]. 
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Fig. 2: plots of EoS parameter u(t) and deceleration parameter, q(t), for F(R) = R 2 , 
f = f e ^(*), / = -0.5 and 6 = 5. Initial values are (f)(0) = -5, 0(0) = 15, ^(0) = 15, 

^(0) = 10, H(0) = 8 and H(0) = 1. 



-0.5-1 1 

(0.28,-1) 




Fig. 3: plot of EoS parameter u eff (z) for F(R) = R 2 ,f = f e b(t>{t \ f Q = -0.5 and 6 = 5. 
Initial values are 0(0) = -5, 0(0) = 15, ^(0) = 15, ip(0) = 10, H(0) = 8 and H(0) = 1. 

4.1. The difference in the distance modulus, fj,(z) 

The difference between the absolute and apparent luminosity of a distance object is given 
by, n(z) = 25 + 51og 10 di,(z) where the Luminosity distance quantity, d L (z) is given by 

d ^ = ^[w)- (16) 

In our model, from numerical computation one obtains H(z) which can be used to evaluate 
n(z). In Fig. 4, we compare fi(z) in our model with the observational data. As can be seen 
the model for the given conditions of dynamical variables is in good agreement with the 
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Fig. 4: The graph of distance modulus fi(z) plotted as function of redshift, 
in our model ( gray curve) in comparison with the observational data (the black data points). 
ICs. 0(0) = -0.5, 0(0) = 0.13, V(0) = -0.6, ^(0) = -0.2, H(0) = 0.1 and H(0) = 1. 

observational data. 



4.2. CPL, CRD and our model 



Following 24J, in Chevallier-Polarski-Linder (CPL) parametrization model one can use 
linearly approximated EoS parameter, 



du 



UJ C pl « W 



cpl 



(a — 1) = u + tor. 



(17) 



where u is current value of the EoS and cui 



da y 1 1 + z ' 

is its running factor. Using the above 



duj, 



cpl 



equation we can find the following equation for Hubble parameter 
H(z) 2 



HI 



n m {\ + z f + (i - n m )(i + z f(^+^ 



x exp 



-3u>i( : 



(18) 



l + z' 

In CPL model the parametrization is fitted for values of c^o, ooi and Q m corresponding to 
our model. On the other hand, the CRD can be extracted from 



:i + z)H -H(z) 



(19) 



which is known as Mc Vittie equation. This equation immediately leads to velocity drift 

cH(z) 



v = cHi 



o 



l + z 



(20) 
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Fig. 5: The graph of velocity drift plotted as function of redshift z for both CPL 
parametrization and our model. ICs. (f)(0) = —0.5, (f>(0) = 0.13, "0(0) = —0.6, 
ip(0) = -0.2, H(0) = 0.1 and H(0) = 1. 

By using equation f lT8|) in CPL model, the velocity drift with respect to the redshift can be 
obtained against observational data. Fig. 5, shows the velocity drift in terms of redshift 
against observational data for our model. 

From Fig. 5, a comparison between our model and CPL model with respect to the 
observational data is given. As can be seen that our model for the redshift z greater than 
3.5 better fits the data than CPL parameter. 

5. SUMMARY 

In this paper, we have studied the evolution of the gravitational fields both analytically 
and numerically in the non-local F(R) model in which the geometry partly modified by two 
scalar fields, one of them a Quintessence-like field, the other a Phantom-like field. 

We discussed a bouncing non-singular cosmology, with an initial contracting phase which 
lasts until to a non- vanishing minimal radius is reached and then smoothly transits into an 
expanding phase which provides a possible solution to the singularity problem of Standard 
Big Bang cosmology, a problem which is not cured by scalar-field driven inflationary models. 
The evolution of EoS parameter, hubble parameter, scale factor, and deceleration parameter 
is numerically obtained. The violations of the null energy condition required to get a bounce 
are obtained for the model, which allows a transition of the EoS parameter through the 
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cosmological constant boundary. In comparison with our previous work in [lj], we find that 
in the analytical discussion of the phantom crossing behavior of the EoS parameter, we need 
also constrain the scalar fields and their first and second derivatives. Moreover, we have 
a,o additiona. constraints on hubb.e parameter and its first and second derivatives. In Q, 
with reconstruction of function /(</>), we found that its behavior is similar to exponential 
function. Here, we are motivated to assume an exponential behavior for f ((/)). Beside, such 
exponential functions have been known lead to interesting physics in accelerated expanding 
cosmological models 251. 

n 

In comparison to our previous work also two cosmological tests are performed to 
validate our model. The distance modulus obtained from numerical computation fitted with 
the observational data with a phantom crossing behavior for the universe is predicted at 
redshift z = 0.29. With regards to the assumed functionality for F(R) and /(</>), and the 
crossing behavior for the model, the second test is performed. From numerical calculation, 
the drift velocity is obtained for the model and is compared with the observational data 
and CPL parametrization. The result shows that our model may better fit the data in 
comparison with CPL model in higher redshifts . 
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